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INTRODUCTION 


Passive  acoustics  are  used  to  detect,  classify,  and  local¬ 
ize  signal  sources.  Present  practice  is  based  primarily  on  spectral 
and  cross-spectral  analysis  in  which  phase  relations  between  signal 
components  of  differing  frequency  is  not  exploited.  Bispectral 
analysis  is  an  important  extension  of  power  spectral  analysis  which 
makes  use  of  inter-frequency  phase  information.  MacDonald  has  been 
promoting  the  application  of  bispectra  in  various  problems  (see  JSR- 
82-601,  Speech  Research,  by  Despain,  MacDonald,  and  Rothaus).  This 
collection  introduces  the  bispectrum,  and  other  polyspectra,  and 
summarizes  three  preliminary  studies  exploring  application  to 
passive  acoustic  ASH. 

The  first  section  is  a  primer  by  G.  MacDonald  on  blspectra 
and  higher  order  spectral  constructs.  In  this  section,  the  auto- 
blspectrum  and  cross-blspectrum  are  Introduced  and  related  to  third 
order  time-lagged  mean  products  of  the  time  series  and  to  higher 
order  transfer  functions  connecting  the  process  under  examination  to 
a  serially  Independent  generating  function.  Also  discussed  is  use 
of  blspectra  in  detecting  and  describing  nonlinearity  in  the  process 
generating  the  time  series. 
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In  Che  second  section,  R.  Davis  addresses  the  detectability 
of  a  signal  through  its  spectrum  and  its  bispectrum.  Surprisingly, 
no  processing  gain  derives  from  bispectral  analysis  even  though  it 
makes  use  of  waveform  (inter-frequency  phase  relation)  informa¬ 
tion.  Reliable  detection  requires  more  signal  energy  for  bispectral 
detection  than  for  spectral  detection  unless  the  signal  skewness  is 
large  compared  with  unity.  The  amount  of  averaging  applied  to 
reduce  random  sampling  variations  is  proportional  to  the  time- 
bandwidth  product  N;  as  N  increases,  the  disadvantage  of  bispectral 
analysis  increases. 

In  the  third  section  H.  Abarbanel  describes  some  experi¬ 
ments  in  which  numerically  generated  signals,  rich  in  harmonic 
content,  were  subjected  to  power-spectral  and  bispectral  analysis. 
These  show  that  the  signal  bispectrum  contains  more  structure  than 
the  power  spectrum.  Indicating  that  bispectral  analysis  of  signals 
may  provide  very  useful  additional  classification  beyond  that 
obtained  from  energy  methods.  Consonant  with  the  results  of 
Section  2,  it  is  found  that  background  noise  obscures  the  signal 
bispectrum  and  that  the  energy  slgnal-to-noise  ratio  required  for 
useful  signal  identification  is  comparable  for  spectral  and 
bispectral  analysis. 


In  the  fourth  section  W.  Hunk  outlines  an  application  in 
which  the  inter-frequency  phase  information  in  a  class  of  signals 
one  might  expect  to  be  produced  by  submarines  can  be  used  to  obtain 
localization  information.  In  addition  to  presenting  this  poten¬ 
tially  practical  application,  this  section  provides  insight  into  the 
origin  and  description  of  inter-frequency  phase  relations  using 
blspectra  and  related  analytic  tools. 

To  summarize  our  preliminary  conclusions: 

(1)  It  is  possible  that  inter-frequency  phase  can  be  used 
to  describe  aspects  of  source-receiver  geometry,  and 
this  may  have  real  utility  in  short  range  target 
trailing. 

(2)  The  bispectrum  magnitude  and  phase  contains  informa¬ 
tion  about  the  signal  not  available  in  the  power 
spectrum.  This  could  provide  useful  classification 
signatures  if  source  waveform  is  not  strongly 
dependent  on  geometry  and  operating  conditions. 

(3)  Blspectra  do  not  appear  to  offer  any  advantage  in 
detectability  at  low  slgnal-to-noise  ratio  unless  the 
signal  skewness  is  large. 


From  a  broader  perspective ,  bispectral  analysis  is  just  one 


analytic  tool  which  describes  inter-frequency  phase  relations  or 
waveform.  These  should,  in  general,  provide  additional  signal 
description  and  some  discrimination  of  signal  and  noise.  To  assess 
this  possibility,  it  is  necessary  to  understand  what  waveform 
signatures  exist  in  real  signals,  how  well  these  survive  propaga¬ 
tion,  and  how  they  may  best  be  analyzed.  Toward  this  end  we 
recommend: 

(1)  Submarine  and  surface  ship  signals  should  be  analyzed 
for  bispectral,  and  other  waveform-descriptive, 
signatures.  Sonobouy  records  obtained  at  high  signal- 
to-noise  ratio  would  be  best  for  this  purpose. 

(2)  Analysis  of  inter-frequency  phase  distortion  during 
propagation  should  be  explored  analytically  to  provide 
order-of-magnitude  estimates. 

Beyond  passive  acoustics,  we  hypothesize  that  bispectral 
analysis,  or  some  other  analysis  tool  based  on  inter-frequency  phase 
relations  may  be  useful  in  active  acoustics  in  the  presence  of 
reverberation.  Assuming  that  the  target  is  relatively  localized 
compared  with  the  reflectors  producing  reverberation,  it  is  likely 
that  phase  information  may  help  discriminate  target  and  clutter. 
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A  PRIMER  ON  HIGHER-ORDER  OR  POLYSPECTRA 


Analysis  of  time  series  is  the  central  problem  of  data 
analysis  in  a  wide  variety  of  defense-  and  intelligence-related 
applications.  Development  of  the  theory  of  time  series  analysis  has 
followed  two  paths,  communications  engineering  and  mathematical 
statistics.  The  engineering  community  has  emphasized  the  frequency 
or  spectral  approach,  while  the  statisticians  have  found  comfort  in 
the  time  domain.  Both  groups  have  based  their  work  on  stationary 
linear  models.  These  models  have  had  great  success  in  a  wide 
variety  of  unrelated  fields,  and  quite  naturally  there  has  been 
little  effort  to  explore  non-linear,  non-stationary  models.  The 
very  great  achievements  and  the  refined  methods  used  in  stationary 
linear  models  suggest  that  it  is  unlikely  that  further  progress  in 
time  series  analysis  is  likely  by  pursuing  concepts  in  areas  such  as 
linear  predictive  coding.  Further,  simple  observations  of  the  world 
around  us  show  that  actual  phenomena  are  non-linear  and  non- 
stationary.  Progress  in  data  processing  is  likely  to  be  in  the 
direction  of  models  which  incorporate  non-stationarity  and  non¬ 
linearity.  This  primer  is  concerned  only  with  non-linearity.  Non- 
stationarity  will  be  dealt  with  in  a  separate  paper. 


Linear  and  Non-Linear  Models 


Any  discrete  stationary  process  xt  with  a  spectral  density 
can  be  represented  as  a  linear  combination  of  an  uncorrelated 
stationary  process  in  the  form 


x 


t 


h 

n 


t-n 


(1.1) 


Equation  (1.1)  provides  a  general  representation  of  a  linear  process 
and  can  easily  be  transformed  into  the  familiar  ARMA  (autoregressive 
moving  average  model)  if  the  function 


g(z)  -l  hn  zn  (1.2) 

n*-°° 

can  be  approximated  by  a  rational  function  a(z)/B(z)  .  Provided 
that  (1.2)  holds,  and  assuming  $(z)  has  constant  term  one,  then 


\  +  6iVi  +  •  •  •  •  Vt  +  Vt-i  +  Vt-2  + 


(1.3) 


in  which  is  an  uncorrelated  process 


E[e  e  I  -  5 
1  t-n  t-mJ  n-m 


(1.4) 


2  2 
where  e  is  the  expected  value  of  e  . 
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This  representation  does  not  distinguish  between  uncorrelated  and 
independent  variables ;  a  linear  model  is  one  in  which 


x 


t 


h 

n 


t-n 


(1.5) 


where  et  are  independent  variables.  Of  Interest  is  that  ,  an 
uncorrelated  process,  and  et,  an  independent  process,  have  the  same 
second-order  structure.  Both  have  white  spectra  but  may  differ  in 
other  ways,  as  will  be  noted  below.  If  e^.  is  a  Gaussian  process, 
it  is  both  uncorrelated  and  independent,  and  the  distinction  between 
and  e£  vanishes;  xt  is  then  also  a  Gaussian  process. 


In  order  to  illustrate  the  difference  between  et  and 
in  general,  we  consider  the  prediction  problem.  If  the  sequence  et 
is  strictly  independent,  then  the  past  contains  no  Information  about 
the  future,  and  the  best  predictor  of  efc  is  its  mean.  This  is  also 
true  for  an  uncorrelated  process  provided  the  predictors  are 

linear .  However,  the  past  may  contain  information  about  the  future, 
which  will  be  revealed  if  the  predictors  are  non-linear  functions  of 
the  observed  values,  even  though  the  values  of  are 

uncorrelated.  For  example,  consider  a  process  n  defined  by 


nt  “ 


et  * 


a  e 


t-1  t-2 
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where,  as  usual,  efc  represents  a  serially  Independent  process  with 

2 

zero  mean  and  constant  variance,  o  .  The  process  n  is 
uncorrelated  and  as  far  as  Its  second-order  properties  are 
concerned,  it  behaves  as  an  independent  process.  However,  unlike 
strictly  Independent  processes,  the  optimum  mean  square  error 
predictor  which  Is  at  most  quadratic,  and  looks  two  steps  back  in 
time.  Is  8 imply 


"t+i  ■  Bnt  Vi 


where  0  -  o/(l  +  a^o^)  . 


From  this  simple  observation,  it  Is  obvious  that  there  are  processes 
that  do  not  follow  the  linear  representation  given  by  (1.5). 


A  generalization  of  (1.5)  is 


L h*  e‘-i + 1 1  hij  et-i  et-j + h«k  et-i  vj  et-k +  •  •  • 

(1.6) 


. 
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in  analogy  Co  a  Taylor  series  expansion.  Such  functional  expansions 
were  first  studied  by  Volterra  and  introduced  into  non-linear  sta¬ 
tistics  by  Weiner,  so  the  expansion  in  (1.6)  is  generally  known  as 
the  Volterra-Weiner  expansion.  One  may  always  suppose  that  the  co¬ 
efficients  hjj  Q  are  unaltered  by  permutation  of  the  indecies. 
Though  there  is  a  growing  literature  on  such  representations,  the 
problem  of  estimating  the  generalized  transfer  functions  hjj  n  has 
proven  to  be  intractable  primarily  because  of  the  large  number  of 
parameters  Involved.  Progress  in  using  such  a  representation  is 
likely  only  if  the  process  can  be  represented  by  a  small  number 
of  parameters  or  if  the  coefficients  hj,j  n  have  some  sort  of 
"smoothness”  property.  In  the  linear  case,  the  smoothness  condition 
is  imposed  on  Hj(f),  the  Fourier  transform  of  h^,  by  insisting  that 
Hj(f)  be  a  decent  function  or,  in  the  case  of  an  ARMA  process,  that 
Hj(f)  arises  from  a  rational  function. 

Since  transfer  functions  have  proven  to  be  exceedingly  use¬ 
ful  in  linear  problems,  an  obvious  generalization  is  to  define  a  set 
of  generalized  transfer  functions  by 


This  is  the  familiar  result  that  a  non-linear  process  acting  on  a 
sinusoidal  input  at  a  frequency  fQ  produces  an  output  containing 
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integer  multiples  of  the  frequency  fQ. 

1.2  Polyspectra 

For  linear  problems,  the  second  order-moments  of  the  exci¬ 
tation  e£  and  the  process  x£  are  used  to  estimate  the  coefficients 
in  an  ARMA  representation,  though  they  do  not  completely  determine 
the  coefficients.  For  non-linear  processes,  higher-order  moments 
and  their  Fourier  transforms  may  provide  an  insight.  For  the 
process  xt,  the  third  and  fourth  cumulants  are  defined  by 


R3^T  1  ,T2 ,T 3 )  * 


R4(T,.t2.t3,,4)  '  *  «.)(*■ 


t-K1  X  '  '  t+tj  X  *  '  t+T  3  X/v  t+T^ 


)  -  O]  E[(^  -  O  -  O 


E^Xt-K  "  V  V"t-K2  "x  ‘ J  “L  v”t+r3  V"t-K4  X 


El^xfKl  -  (Xt-K3  “  E  lK-K2  -  (xt+r4  - 


Ei(*^  -  O  -  Oi  E  lOw  -  O  (x—  - ' 


•t+Tj  ~x'  V"t+I4  -x'M  ~  t+c2  x'  v  t+r3 


*  •  *> 


where  Is  the  nean  of  xt 


This  notation,  while  unconventional,  clearly  reveals  the  syn try  in 
the  t 's  .  We  note  that  for  a  stationary  process,  the  origin  can  be 
selected  so  that  a  given  t  can  be  set  to  aero,  e.g., 
in  which  case  (1.9)  reverts  to  the  usual  fora. 


The  third-  and  fourth-order  polyspectra  are  then  defined  by 


00  00 


Bj(fj»f2»fj)  “X  I  I  * 


-2*i  (f  j+f  2t  2+f  3T  3  ) 


•  t2«- “  Tj 


(1.10) 
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The  Inverse  transforms  are  then 


R3(VVT3) 


V2  v2  v2 

-V2  -V2  -v2 


2iri(f1T1+f2t2+f3T3) 
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and  similarly  for  R4(T1»T2»T3»T4)  • 


The  relationships  among  the  frequencies  in  B3  and  follow 
from  the  definition  of  and  R4.  For  the  representation 


'2  2wif  t 
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then  R^  becomes 
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For  historical  reasons,  B  Is  known  as  the  bispectrum  and  the 
trlspectrua  even  though  (1.13)  suggest  that  B  should  really  be 
called  the  trispectrum  and  so  on;  we  will  abide  by  present  conven¬ 
tion. 

1.3  Polyspectra  for  Gaussian  Processes 

A  well-known  result  of  statistical  theory  is  that  all  joint 
cumulants  higher  than  the  second  order  vanish  for  a  multivariate 
normal  process.  The  immediate  and  Important  consequence  is  that  all 
spectra  of  order  higher  than  two  vanish  for  a  Gaussian  process.  The 
higher-order  spectra  would  thus  appear  to  be  of  use  In  investigating 
properties  of  a  non-Gaussian  stochastic  process  or  of  a  non-linear 
system  driven  by  a  random  input.  As  noted  above,  polyspectra  give  a 
measure  of  the  phase  correlation  between  components  whose  frequen¬ 
cies  sum  to  zero. 

While  the  above  paragraph  captures  the  essential  features 
of  polyspectra,  there  have  been  few  genuine  applications.  In  part, 
this  has  arisen  from  the  large  data  volumes  and  heavy  computation 
required  to  numerically  evaluate  higher-order  spectra  and  in  part 
from  statistical  difficulties  in  interpreting  the  results.  However, 
it  would  appear  that  the  most  important  impediment  has  been  the  lack 
of  a  physically  understandable  interpretation  of  polyspectra.  As  a 
first  step,  we  consider  cross  polyspectrum. 
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The  cross-bispectrum  is  a  simple  generalization  of  the 
bispectrum.  A  form  of  the  third-order  moment  for  two  processes  y^ 
and  xfc,  where  for  tutorial  purposes  xt  can  be  thought  of  as  an  input 
into  a  system  and  yt  the  output,  is 


Wvv'j)  ■  'to 


t+T 


and  the  corresponding  blspectrum  is 


Bxxy^fl’f2*f3^  \L.L.  L.  Rxxy^Tl’T2’T3^  e 


-2iri(f1x  1+f2T2+f3T3 


x^m-co  T^— 


and  if 


v2 


2xift 

i  < 

-1/,  7 


/  .  e‘-"“‘-dZ,.(f) 


then 


B«y<fl-£2'f3>d£ldf2df3  -  2[dZx(f1)d2!t(£2)d2y(f3)] 

An  illustrative  example  of  the  use  and  limitations  of  the 


bispectral  analysis  is  contained  in  the  process 


(1.16 


-  I  hlXt-£  +  l  I 

1-0  1-0  m-0 


h  x  .x  +  N 
In  t-1  t-m  t 


where  Is  a  zero  mean  stationary  Gaussian  process,  and  is  a 
zero  mean  noise  independent  of  xfc.  From  (1.8),  we  see  (1.16)  takes 
the  form 


yt 


/  ^  H  (f)  e2lliftdZ  (f) 

-v2  1 


(i.i7: 


V2  V2  2iri(f.+fjt 

+ ;  i/ ;  i/  H2(fi*f2]  e  dMfiK(f2) +  Nt 

-v2  -v2 


The  mean  of  yt  is 


/  H-(f ,-f )P(f )df 

-v2  2 


(1.18 


where  P  is  the  power  spectrum 

P(f)df  -  E[dZ  (f)dZ  (-f)] 

The  first  odd-order  moment  of  yt  is  thus  dependent  on  the  quadratic 
transfer  function  H-(f.,f-)  . 


The  second -order  moment  E[y  yt+T] 


takes  the  form 


VT)  *  Etv  yt«l 

V2  V2  2«i(f.t+f  t )  , 

-  /  *  /  l  ■,(*,)>,(*,)  •  1  2  '2,lfI  yK(^K('2)] 

“  '2  “72 


V2  V2  V2  V2  2iri(f  +f  )t 

+  /  /  /  /  H2(f 1«f2^  H2(f3,£4)  e 

“72  “V2  -  v2  -v2 


2*i(f  +f  Jt 

*  e  E[dZx(f1)dZx(f2)dZx(f3)dZx(f4)] 

(since  x  is  Gaussian,  E[dZ  (f.)dZ  (f  )dZ  (f_)]  *  0  ).  Using  the 
standard  decomposition  theorem  for  even-order  moments  of  Gaussian 
processes , 


E(dZ*(flK(f2K(f3>dZ*(£4)]  * 


E [dZx (f  j }dZx (f  2 ) ]E [dZx (f 3 )dZx (f 4 ) ] 

+  E[dZx (f  L )dZx (f 3 ) ]E [dZx (f 2 )dZx (x4 ) ] 


+  E(dzx  (f  !  )dzx(f  4 )  ]E  [«x(f  2  )dzx(f  3 )  1 


and  (1.17)  reduces  to 


Vo 

R  (t)  -  /  H.(f)H  (-f)  e2irifTP(f )df 

yy  -Vo  1  1 


+  /  J  H2(flf-f 1)H2(f3,-f3)p(fl)p(f3)df tdf3 

-  y2  - 1/2 


Vo  Vo  -2iri(f  x+f  t) 

+  /  ,  /  :  e  P(f1)p(f2)df1df2 

-Hi  ~  V2 


Vo  Vo  -2iri(f  T+f.x) 

+/  :  /  .  H2(f1,f2)H2(-f2.-f1)  e  p(f1)p(f2)df1df2 

-V2  -v2 

The  second-order  moment  of  yt  thus  contains  products  of  the  power 
spectra  of  x^.  weighted  by  the  quadratic  transfer  function. 

An  examination  of  the  cross  bispectra  is  more  revealing  of 
the  underlying  structure  of  the  generating  process.  The  cross 
moment  is 


Ry*(n)  "  E[,t+nKtl  ’  E!^0  hixtxt-frtu  +  J^g  hSmxt-i-Hixt-iB-HixC ^ 


I  ‘  .1  htRxx(n'*) 


since  the  odd-order  moments  of  the  Gaussian  process  vanish.  Taking 
the  Fourier  transform  of  both  sides,  we  obtain 


B^Cf)  -  H1(f)Bxx(f)  -  H1(f)P(f) 


and  arrive  at  the  familiar  result  for  linear  systems: 


Vf)  "  Byx(f)/P(f) 


The  cross-bispectrum  B _  can  be  obtained  from  the  representation  of 

yxx 

yt  given  in  (1.17): 

f  ,  v2  ,V2  V2  2U[£1t+f2(t«)+f3(fHi)] 

•  H,(f)E[<lZ]i(f1)dZx(f2)d2i(f3)j 

,V2  V2  V2  V2  2.l[(f  «  )t«  (t+l)+f4(tt.)] 

+  /,/,/,/,«  H  (f  *  ) 

-v2  -v2  -v2  -v2 


E(dZx(f1)dZit(f2)dZxff,)dZx(f4)] 


The  first  Integral  vanishes  since  x£  Is  Gaussian,  and  the  second 
tern  can  be  decomposed  as  above  so  that 


'tVt+.VJ  -  I  t  H2(£l'-fl)P,(fl)Px(f3)  « 

-  V2  -  V2 


2rif3(l-m) 


df 


l/2  l/2  2iri(lf  ,+mf-) 

+  /  „  /  E2(-£2--fl>x(£l)Px(f2)  ‘ 

-  12  ~  '2 


dfldf2 


V2  V2  2iri(f  m) 

+  /  „  /  l7  H2(-fl--f2)Px(£l)P,(f2)  * 

-v2  -v2 


dfldf2 


We  now  consider  the  quantity 


Vx<1"*>  '  “frtWwJ  ' 

From  (1.18)  we  note  that  the  first  integral  In  (1.20)  is  just 

.  •«  ch.t 

V2  V2  2iri  ( If .  +mf  y ) 

Ryxx(‘-*>  ■  2  /  „  /  „  H2(-£l,£2)p(£1)p(£2)  .  «,« 


Since  h 


can  be  made  symmetric  in  t  and  m,  then  ^(fj.fj) 
symmetric  in  fj  and  f2*  The  cross-bispectrum  Byxx(fi»f2^  is  then 

and  the  quadratic  transfer  function  is  determined  by 

■  V*K--£2)/2P(fl>P(£2> 

For  the  particular  process  given  in  (1.17),  it  is  possible  to 
recover  the  linear  and  quadratic  transfer  functions  by  measuring  the 
power  spectrum  of  the  input  xt>  the  cross  spectrum  ,  and  the 
cross-bispectrum  ®yXX(^>^2^  *  ^or  ®°re  Senera^  processes  with 
terms  of  various  orders,  this  is  no  longer  possible,  and  an  itera¬ 
tive  calculation  is  required. 

1.5  Meaning  of  a  Finite  Bispectrum 

If  a  time  series  yt  has  a  statistically  significant  non¬ 
vanishing  bispectrum,  then  two  Interpretations  are  possible.  The 
generating  process  contains  non-linearities  which  may  be  described 
in  the  form  of  (1.6)  or  possibly  some  other  representation,  with  the 


excitation  function  being  Gaussian  or  non-Gaussian .  Alterna¬ 
tively,  the  generating  model  for  the  process  is  purely  linear  but 
the  excitations  efc  are  non-normal,  have  a  finite  skewness,  and  are 
independent.  In  this  case,  the  third-order  moment  for 


‘t  ■  i.Vt« 


*  Et6tl  l  htht+»ht-H. 


since  efc  are  independent.  From  the  definition  of  the  linear 
transfer  function  Hj(f)  , 


H  (f)  -  l  he 
t-0 


i2irft 


the  bispectrum  of  xt  is 


Bx(fi*f2)  “  Elet^  Hl(fJ  M£2J  '  f2^ 


The  power  spectrum  of  xfc  is  just 


P^(f)  -  const  Hj(f) 


so  that  the  square  of  the  bicoherence  given  by 

»*(*,)  *x(f2)  px(fi  +  f2> 

is  also  a  constant.  Thus  for  a  linear  system  with  a  Gaussian  input, 
the  bicoherence  will  vanish,  and  a  linear  system  with  a  non-Gaussian 
input  will  have  a  constant  bicoherence  not  zero,  provided  the  input 
has  a  finite  skewness. 

The  determination  of  the  bispectrum  does  not  by  itself 
permit  the  identification  of  a  non-linear  underlying  process  or  of  a 
non-Gaussian  excitation  function.  The  excitation  function  could  be 
non-Gaussian  but  have  a  vanishing  skewness,  and  the  resulting 
blspectra  of  xt  would  vanish.  Alternatively,  a  constant  bispectrum 
does  not  imply  a  linear  process.  A  finite  and  varying  bispectrum  is 
certainly  suggestive  of  non-linearities  in  the  generating  process. 

It  is  for  this  reason  that  the  bispectra  analysis  of  speech,  with 
strongly  bispectral  peaks  for  the  vowels,  suggests  that  speech 
production  has  Important  non-linear  characteristics.  Similarly,  in 
underwater  acoustics  bispectral  signatures  may  provide  powerful 
methods  for  classifying  sound  sources. 
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BISPECTRAL  DETECTION 


Introduction 


Present  practice  in  acoustic  detection  rests  on  dlscrlmi- 
nating  signal  from  noise  on  the  basis  of  energy  distribution  in 
frequency  and  time.  That  is,  waveform  and  phase  structure  of  sig¬ 
nals  are  not  exploited.  The  existence  of  sequences  of  harmonically 
related  lines  in  signal  spectra  suggests  that,  at  least  in  the  case 
of  machinery  produced  sound,  there  may  be  useful  phase  information 
in  signals  of  Interest.  Beyond  this,  it  might  be  supposed  that 
broadband  signals  produced  by  highly  nonlinear  hydrodynamics  (flow 
noise  and  blade-rate  signatures)  are  also  associated  with  stable 
phase  relations  between  signal  components  of  differing  frequencies. 
The  question  addressed  here  is  whether  such  phase  relations,  if  they 
exist  and  are  not  destroyed  by  propagation,  could  be  used  to  enhance 
detectability  of  weak  signals. 


He  consider  the  received  record,  R(t),  to  be  the  sum  of  a 
signal,  S(t),  and  noise,  U(t).  Over  any  record  of  length  T  starting 
at  time  t  these  have  Fourier  representations 


[S(t);  U(t)]  -  I  (s(w,t);  u(w,t)J  e 
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where  u  is  an  integer  aultiple  of  Aw  «  2*/T  and 


[s(w,t ) ;  u(w,x )]  /  (S( t) ;  0(t)]  e~iU,tdt  . 

The  received  record  has  Fourier  amplitudes 

r(w,f)  -  s(w,t)  +  u(ui,t ) 

For  this  discussion  it  will  be  assumed  that  records  are 
prefiltered  so  that  the  power  spectrum  of  noise  is  uniform,  that  is 

pu<“>  =  li“Au>-*o<lu<“,T)|2>/A“  *  Pu 

where  <  >  denotes  a  long-time  average  over  T  .  Further,  it  is 
assumed  that  noise  is  the  sum  of  many  contributions  from  independent 
sources  and  may  be  approximated  as  normally  distributed.  Neither  of 
these  idealizations  is  strictly  accurate  for  oceanic  acoustic  noise, 
but  they  are  fair  approximations  which  do  not  prejudice  the  compari¬ 
son  made  here. 

The  interest  here  is  in  learning  how  use  of  additional 
information  about  the  nature  of  the  signal,  S(t),  affects  detect¬ 
ability.  For  this  comparison  three  simple  detectors  are  considered. 


Energy  detection  Is  based  on  some  variant  of  the  spectro- 


graa 


P(w,t) 


where  the  average  received  power  Is 

t+T 

Y  /  R  (t)  dt  -  i  is  p(«,t) 

X  (I) 

The  spectrum  is  the  time  average  P(oi)  -  <P(w,x)>  .  The  simplest 
energy  detector  is  of  the  form 

a  °u 

D_(x)  -  Z  Au>  P((i),T )  -  ~  Z  Au  P(u>,x)  (2. IE) 

c  n 

u>  u  u> 

where  the  sums  are  over  particular  frequency  regions  ft  and  fl^  . 

Here  ft  is  a  region  of  positive  frequencies  occupied  by  the  signal 

and  ft^  is  a  non-overlapping  region  in  which  signal  is  absent.  N 

and  Nu  are  the  number  of  fundamental  frequency  intervals  in  the 

respective  frequency  regions;  N  is  the  time -bandwidth  product 

ZAo>  •  T/2ir  .  The  frequency  regions  need  not  be  continuous  and  ft 
ft 

might,  for  example,  be  a  series  of  bands  encompassing  various 


anticipated  Machinery  signal  lines.  The  second  term  in  (2. IE) 


siaply  serves  to  subtract  froa  D  an  estlaate  of  the  noise 

IS 

contribution  in  the  first  tera.  In  considering  detector  output 
variability,  it  is  assuaed  that  Ny  »  N  . 


Blspectral  detection  might  be  based  on  the  bi-spectrogram 

_  r(u  ,t)  r(w  ,t)  r(-aj  -  u>  ,t) 

^ - L-J— 

(Ato)2 

where  the  received  signal's  average  cube  is 

t+T 

Y  j  R3(t)dt  -  Z  Aw  E  Aw  BCWpWj.T)  . 

T  (i)  hs 

1  2 

The  bispectrum  is  the  time  average 

_  !♦(«..  O 

B(o»ltw2)  -  <B(u)lf(o2,T)>  -  |B(a1,«2)|  e  , 


where  $  is  the  bi -phase.  From  the  view  of  signal-noise  discri¬ 
mination,  the  bispectrum  is  a  simple  construct  which  makes  use  of 
the  phase  differences  between  different  frequencies;  the  particular 
phase  information  preserved  is  that  which  contributes  to  skewness. 


A  simple  bispectral  detector  is 

A 

8  8  ,  „  -i<Ko>,  .O 

D  (t)  -  Z  Z  (Awj  B(u>  ,(D.  ,t)  e  (2. IB) 

-1-2 

A 

where  $  is  an  a  priori  estimate  of  the  bi -phase  $  .  Note  that 
since  includes  the  Fourier  amplitude  at  frequency 

o>j+  o>2  ,  this  detector  makes  use  of  energy  outside  8  .  Use  of 
this  detector  requires  that  the  bi-phase  be  somewhat  predictable  so 
that  the  various  terms  in  the  sum  (IB)  can  be  phased  to  interfere 
constructively.  In  practice,  detection  could  be  based  on  |Dg|,  in 
which  case  the  only  requirement  would  be  that  be  reason¬ 
ably  constant  for  ui^  and  In  8  . 

Maximal  signal  waveform  Information  utilization  is  embodied 
in  a  detector  based  on  matched  filtering  with  a  perfect  replicate  of 
the  signal.  Such  a  detector  can  be  expressed  as  a  convolution  of 
signal  and  the  record  and  is  equivalent  to 

a 

DM<t)  -  E  8(u>,t)  r(-ai,T)  (2.1M) 

m 

This  is  not,  of  course,  a  realizable  detection  option,  since  S  is 
not  known;  it  is  included  here  only  for  comparison  purposes  to  show 
what  could  be  achieved  if  all  signal  waveform  information  could  be 
used. 
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The  three  detectors  have  coaon  features.  All  involve  a 
sample  of  information  from  a  region  of  frequency  space  (2  and  time 
T;  if  Aft  is  the  bandwidth  of  ft  ,  the  size  of  this  sample  is  th 
time -bandwidth  product  N  -  Aft  •  T/2*  .  If  applied  to  a  record 
consisting  of  only  noise,  all  detectors  will  produce  some  output 
D'(t)  but  the  average  output  over  many  realizations, 

<D'(x)>  ,  will  vanish.  In  the  case  of  energy  detection, 

<DI>  -  0  because  a  priori  information  about  the  noise  spectrum  was 

Ei  “*— 

used  to  make  <D'>  «  N  P  -  rr  N  P  -  0  ;  in  practice  this  may  not 
E  u  N  u  u 

u 

be  strictly  achievable  but  the  bias,  <D^>  ,  will  be  small. 

<D'>  *  0  because  the  bispectrum  of  normally  distributed  noise, 

B 

having  no  stable  inter-frequency  phase  relations,  vanishes;  in 
reality  the  bias  may  not  strictly  vanish  for  some  noise  sources  but 
it  should  be  small  unless  the  noise  arises  from  a  localized  source 
which  we  here  consider  a  signal.  <D*>  *  0  because  noise  is  not 
correlated  with  signal. 

2.3  Simplified  Performance  Measures 

Detector  sensitivity  to  signal  may  be  characterized  by  the 
mean  output,  <D>.  For  comparison  purposes,  the  value  of  this  mean 
output  is 
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o 

where  E  «  <S  >  is  the  signal  variance. 


The  bispectral  detector's  signal  output  depends  on 


ft  ft 


l  l  <s(u1)  s(o>2)  s*(u1  +  t»2)  > 


CO  CO 
1  2 


[  £  <|s(o>)|2  >  ]3/2 


In  the  simplest  case  when  the  bi-phase  $  and  its  estimate  <|>  both 
vanish,  q  is  simply  the  signal  skewness  <S^>/E^^.  When  the  inter¬ 
frequency  relationship  is  more  complex  there  is  no  simple  interpre¬ 
tation  of  q  but  it  remains  true  that  there  is  no  a  priori  bound 
which  can  be  placed  on  its  magnitude.  If  the  probability  density 
function  of  S  is  sufficiently  dispersed  (i.e.  the  probability  of 
extreme  values  is  much  greater  than  for  a  Gaussian  distribution  with 
the  same  variance)  then  large  values  of  q  are  possible.  On  the 
other  hand,  if  the  probability  density  is  central  (with  extreme 


S  V 


values  having  probabilities  comparable  to  a  normal  distribution) 


then  q  will  be  of  order  unity.  Without  knowledge  of  the  signal's 


bi-spectral  character  no  finer  bound  can  be  placed  on  q.  Our 


approach  here  is  aimed  at  the  order-unity-q  case.  We  note  that  if 


the  distribution  of  S  is  dispersed  on  q  can  be  large,  but  to  make 


use  of  the  extreme  values  of  S  which  make  q  large  will  require  use 


of  long  range  records ,  much  longer  than  the  T  required  to  achieve 


frequency  resolution.  Further,  the  extreme  events  responsible  for 


large  q  would  also  be  large  individual  events  to  an  energy  detector 


(such  as  Dg)  and  might  be  detected  on  this  basis. 


In  principle,  (2.2E)  is  valid  for  any  0  ,  T  but  in  prac¬ 


tice  time  varying  sources  or  Doppler  shifting  would  require  ft  to 


increase  with  T.  Achieving  (2.2B)  requires  that  the  phase  of  the 


bl-spectrogram  be  perfectly  stable  and  predicted  by  of  (2. IB). 


Frequency  dependent  propagation  delays  would  thus  degrade  this 


output  in  practice  and,  to  the  extent  that  differential  propagation 


delays  Increase  with  frequency  separation,  would  place  an  upper 


limit  on  the  bandwidth  of  ft  ;  similarly,  time  varying  source  struc¬ 


ture  or  Doppler  shifting  would  make  prediction  of  variations  of  bi¬ 


phase  more  difficult  as  T  increases.  Achieving  (2.2M)  requires  a 


perfect  matched  filter  in  the  presence  of  propagation  phase  distor¬ 


tion  and  the  requirements  for  this  become  more  stringent  as  T  or  the 
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bandwidth  of  ft  increases.  In  the  spirit  of  a  simple  comparison  we 
have  assumed  that  T  and  ft  are  chosen  such  that  frequency  and  phase 
stability  do  not  degrade  the  detector  signal  outputs  (2.2). 


In  the  absence  of  signal  the  detector  outputs,  D',  are 

2  1/2 

random  variables  with  standard  deviations,  o  ■  <|d'|*>  ,  which 

can  be  computed  under  the  assumption  that  the  noise  is  normally 
distributed  and  stationary: 


1/2 

oD  -  H  7  P  Aw 
E  u 


(2.3E) 


Ofi  *  /2  N  [PuAa»] 


3/2 


(2.3B) 


°M  '  El/2lfu4“|1/2 


(2.3M) 


Here  N  is  the  number  of  (positive)  fundamental  frequencies  in  ft  . 
The  probability  density  function  of  each  D'  differs  and  depends  on 
N.  For  N  -  1  (extreme  narrow  band  processing),  D'  is  distributed 

Ci 

2  2  2 

as  X2^X2>  “  1  where  X2  *8  a  chi-squared  variable  with  two 

degrees  of  freedom.  is  distributed  as  Yj^*^  exp(i8)  where 

2 

and  Y2  are  independent  X2  variables  and  0  is  uniformly 
distributed  over  [ 0 , 2ir ]  .  0^  is  normally  distributed,  as  are 
D'  and  D’  in  the  limit  N  «  . 

£  D 
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Comparison  of  the  different  detectors  is  made  difficult 


because  (1)  <Dg>  depends  on  the  slcewness,  q,  and  the  signal  variance 
E  whereas  <Dg>  and  <D^>  depend  only  on  E,  and  (2)  the  bispectral 
detector  makes  use  of  energy  outside  the  signal  bandwith  £2  whereas 
the  other  detectors  depend  only  on  the  energy  within  ft.  There  is 
no  limit  to  the  size  of  the  skewness  and  the  larger  q,  the  more 
useful  is  bispectral  detection. 

For  comparative  purposes,  results  are  presented  in  terms  of 
the  signal  to  noise  ratio 


a 

P  -  E/E  Pu(<d)  Aid  , 


the  ratio  of  signal  variance,  E,  to  the  noise  variance  in 
ft  ,  NP^  Aid  .  In  terms  of  p  ,  the  mean  signal  to  noise  ratios  are 


<V/0E 


(2.4E) 


<DB>/0B  "  <N/2>1/2  <1  P3/2 


(2.4B) 


v  /  — 


..1/2  1/2 


/  A  / 


2.4  Comparison 

To  compare  detectabilities,  we  imagine  the  detection 

criteria 


de  *  dec  ’  8641  <V  >  dbc*  8641  >  dmc  ’ 

where  the  thresholds  DgC,  Dg^ ,  are  set  to  produce  a  specified 
probability  of  false  alarm  (PFA).  The  signal  strength  required  for 
reliable  detection  at  the  specified  PFA  is  then  characterized  by 
p(PFA)  ,  the  signal-to-noise  ratio  required  to  produce  signal 
output  <D>  equal  to  the  threshold  for  false  alarm  rate  PFA. 

Of  fundamental  interest  is  how  the  detectable  signal-to- 
nolse  ratio  depends  on  the  time-bandwidth  product  N  as  N  ♦  •  .  li 
this  case  all  detectors  have  normally  distributed  outputs  and 


Pg(PFA)  -  if 1/2  G(PFA) 

(2.5E 

p  (PFA)  -  q'2/3  if 1/3  G(PFA)2/3 

D 

(2.5B 

p.,(PFA)  -  |  N"1  G(PFA)2 


(2.5M 


where  G(PFA)  is  the  value  a  normally  distributed  variable  with  aero 
mean  and  unit  variance  exceeds  with  probability  PFA.  From  the 
dependences  of  p  on  N  in  (2.5),  it  is  clear  that  the  bispectral 
detector  is  inferior  at  large  N  compared  with  the  energy  detector 
(which  makes  no  use  of  phase  information  in  the  signal)  and  the 
matched  filter  detector  (which  makes  use  of  all  information  about 
signal  waveform).  The  power  of  G(PFA)  appearing  in  (2.5)  depends  on 
which  moment  of  the  noise  contributes  to  the  detector  fluctuation; 

is  essentially  a  first  moment  detector,  Dg  second  moment,  and  Dg 
third  moment.  Since  G(PFA)  is  such  a  weak  function  of  PFA,  the 
power  of  G  in  (2.5)  is  not  a  significant  practical  consideration 
except  at  extremely  small  PFA. 

For  completeness,  the  extreme  of  narrowband  processing  with 
N  ■  1  is  associated  with  the  following  values  of  p  : 


The  values  of  p£  and  are  determined  from  the  X2  an^  normal 

distributions;  p  was  found  by  Monte  Carlo  simulation  using  a 

0 

random  number  generator. 

It  is,  of  course,  possible  to  further  reduce  the  required 
slgnal-to-noise  ratio  for  reliable  detection  by  incoherently 
averaging  the  energy  and  blspectral  detectors.  Use  of  a  detector  of 
the  form 

1  M 

D(m)  -  £  Z  D(t  ) 

M  m-1  “ 

—  1/2 

reduces  the  standard  deviations,  0  ,  in  (3)  by  a  factor  of  M  '  . 

In  this  case,  the  relations  (2.4)  and  (2.5)  pertain  so  long  as  N  is 
the  total  time-bandwidth  product  Mfi/Au  . 

2.5  Conclusion 

The  energy  slgnal-to-noise  ratio  required  for  signal  detec¬ 
tion  by  spectral  and  blspectral  detectors  depends  on  both  the 
detector  time-bandwidth  product,  N,  and  the  acceptable  probability 
of  false  alarm,  PFA.  As  N  increases,  blspectral  detection  shows  a 
relative  disadvantage  of  the  order  Its  disadvantage  decreases 

as  the  PFA  is  reduced.  However,  for  realistic  time -bandwidth 


products  the  bispectral  detection  disadvantage  is  modest.  Even  for 
a  bandwidth  of  100  Hz  and  a  coherent  processing  time  of  2  hours,  the 
factor  corresponds  only  to  a  10  db  difference  in  required 
signal-to-noise  at  a  PFA  of  0.001;  one  minute  Integration  corre¬ 
sponds  to  6  db.  The  signal-to-noise  required  for  bispectral  detec¬ 
tion  also  depends  on  the  signal  skewness,  q.  If  q  is  large  enough, 
it  can  compensate  for  the  bispectral  disadvantage.  Bispectral 

detectability  also  depends  on  how  well  the  theoretical  maximum 
(2.2B)  can  be  approached.  The  critical  optimistic  assumption  in 
(2.2B)  is  that  the  bi-phase  of  the  received  signal  is  predictable. 
This  requires  (1)  that  the  signal  bl-phase  at  the  source  be  stable, 
(2)  that  it  be  predictable,  and  (3)  that  inter-frequency  phase  rela¬ 
tions  not  be  significantly  altered  by  propagation  from  source  to 
receiver. 

We  cannot  assess  the  degree  of  optimism  associated  with  the 
assumption  of  predictable  signal  bi-phase.  First,  we  have  no  infor¬ 
mation  on  the  inter-frequency  phase  relations  in  real  signals.  It 
is  expected  that  there  will  be  stable  phase  relations  between 
machinery  harmonics,  but  if  these  depend  on  geometry  (as  supposed  in 
Hunk's  localisation  concept),  or  on  source  identity  (which  would 
permit  bispectral  classification),  they  will  be  difficult  to  predict 
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3.0 


BISPECTRAL  EXPERIMENTS 


3.1  Introduction 


We  are  interested  here  in  exploring  (1)  the  bispectral 
signatures  of  signals  generated  by  highly  nonlinear,  quasi-perlodic 
processes,  and  (2)  the  efficacy  of  bispectral  analysis  in  extracting 
such  signals  from  noise.  Toward  this  end,  signals  rich  in  harmonic 
content  were  generated  numerically,  purposely  contaminated  by 
Gaussian  white  noise,  and  subjected  to  spectral  and  bispectral 
analysis . 


The  signal  a(k)  k  -  1,  .  .  .  N  was  generated  by  solving  the 
coupled  nonlinear  differential  equations 

da  . 

.  a  -  y  -  sin  z 

-  -y  +  aw  +  ga 


in  which  we  have  one  adjustable  parameter:  g.  As  g  varies,  so  does 
the  topological  character  of  the  solution  to  (3.1).  For  g  <  1.5 
the  solution  for  long  times,  t,  isa*y»w-z*0.  For  g  >  1.5 
periodic  and  more  complicated  asymptotic  motion  ensues. 

We  solved  (3.1)  for  a  variety  of  g  in  the  range  2  <  g  <  4 
for  the  initial  conditions  a(0)  -  y(0)  -  w(0)  ■  1.0  and  z(0)  ■  0.5, 
which  has  no  special  significance.  We  then  stepped  the  solution 
for,  typically,  10^  steps  and  discarded  the  first  few  thousand 
points  to  remove  the  Influence  of  transients  reflecting  the  initial 
conditions.  In  Figure  3.1  we  show  a(k)  -  a(kAt),  where  At  is  a 
fixed  time  step  used  in  the  solution  of  (1),  at  g  -  2.5615.  The 
initial  transient  is  visible  as  is  the  final  asymptotic  behavior. 

The  latter  is  shown  again  in  Figure  3.2  to  emphasize  the  non- 
sinusoidal  nature  of  our  signal.  In  Figure  3.2  a(lc)  minus  its 
average  (av)  is  displayed.  All  results  to  be  shown  here  correspond 
to  the  choice  g  *  2.5615  shown  in  Figure  3.2. 

3.2  Power  Spectrum 

The  power  spectrum  of  the  signal  a(k)  -  av  was  generated 


from  a  series  of  2048  points  taken  from  the  asymptotic  orbit  (Figure 
3.2).  A  cosine  taper  was  applied  to  the  first  fifty  and  last  fifty 


points  to  smooth  the  discontinuities  of  the  chopped  signal.  Figure 
3.3  shows  the  natural  log  of  the  power  spectrum  for  the  first  64 
frequencies  of  the  1024  computed.  In  Figure  3.4  the  averaged  power 
spectrum,  generated  by  taking  2^  -  64  samples  of  length  2**  ■  2048 
from  a  data  sample  of  2^  points  of  a(k),  is  presented.  Each  sample 
was  tapered  and  Fourier  transformed,  then  the  spectrum  was  averaged 
over  64  samples.  This  averaging  should  have  little  effect  on  the 
signal  spectrum,  which  is,  in  fact,  what  we  see.  Averaging  elimi¬ 
nates  some  of  the  effects  of  noise  in  our  samples  (noise  arising 
from  machine  round  off  errors)  and  of  the  phase  relation  between  the 
signal  oscillations  and  the  record  ends. 

Next  we  added  Gaussian  random  noise  with  a  white  spectrum 
to  the  signal  and  processed  the  noisy  signal  in  the  same  fashion. 

In  Figures  3.5  and  3.6  the  power  spectrum  of  the  contaminated  signal 
is  shown  for  two  different  levels  of  power  in  the  noise.  In  Figure 
3.5  everyone  would  see  the  fundamental  and  sharp  eyed  optimists  also 
can  make  out  a  harmonic  (or  two?).  In  Figure  3.6  even  the  optimists 
are  restricted  to  being  fundamentalists. 

The  problem  of  detecting  the  signal  spectrum  peaks  in 
Figures  3.5  and  3.6  is,  of  course,  an  example  of  energy  detection  as 
discussed  in  Section  2.  Within  this  context,  the  task  is  to 
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Figure  3.3  Natural  logarithm  of  the  signal  power-periodogram  computed  from  a  single 
realization.  Only  the  lowest  64  frequencies  are  shown. 
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discriminate  between  spectral  structure  caused  by  the  signal  and 
random  fluctuations  of  the  noise  spectrum.  The  standard  deviation 
of  the  random  fluctuations  given  by  (2.3E)  compares  well  with  the 
evidently  random  wiggles  in  Figures  3.5  and  3.6. 

3.3  Bispectra 

Now  we  turn  to  the  bispectrum.  First  examined  was  the 
diagonal  bispectrum 

B(f,f)  -  <I(f)  a(f )  I*(2f)> 

where  a(f)  is  the  Fourier  transform  of  a(k)  -  av,  with  the  usual 
cosine  taper.  In  Figure  3.7  we  display  log  |B(f,f)|  from  one 
sample  of  2**  points;  64  frequencies  are  shown  and  there  is  no  noise 
added  to  the  signal.  In  Figure  3.8  we  have  taken  64  samples  of  the 
noise  free  signal  and  averaged  B(f,f)  over  that  ensemble.  Since 
B(f,f)  is  cubic  in  the  signal,  averaging  should  reduce  the  effect  of 
any  Gaussian  noise  (round  off  error)  and  there  is  some  evidence  of 
this  in  Figure  3.8. 
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Figure  3.8  Natural  logarithm  of  the  squared  magnitude  of  the  signal  diagonal  bispectrum 
taken  from  64  realizations.  Dominant  peaks  occur  at  f0,  3f0,  5f0 . . .  but  additional 
structure  is  evident  at  1/2  f0,  5fo/2,  9fo/2,  llfo/2. 
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Figures  3.9  and  3.10  show  the  diagonal  bispectrum  of  signal 
plus  Gaussian  noise  for  the  same  noise  levels  as  used  in  Figures  3.5 
and  3.6  for  the  power  spectrum.  It  is  quite  clear  that  an  ensemble 
of  64  members  has  not  produced  sufficient  reduction  of  background 
noise  to  make  the  bispectrum  a  significantly  more  efficient  signal 
identifier  than  is  the  power  spectrum.  This  is  consonant  with  the 
conclusions  of  the  detectability  assessment  in  Section  2.  In  fact, 
the  bispectrum  standard  deviation  of  (2.3B)  corresponds  to  noise 
levels  of  31.5  and  35.5,  respectively,  in  Figures  3.5  and  3.6.  This 
Indicates  that  most  of  the  structure  in  these  figures  is  due  to 
random  sampling  errors. 

We  made  one  further  exploration  into  the  structure  of 
B(f j ,f 2>  •  ®y  locating  the  fundamental,  fQ  ,  of  the  signal  from 
the  power  spectrum  we  calculated 

B(fo,f)  -  <a(fQ)  a(f )  r\i()+f)> 

averaged  over  a  few  frequency  bins  on  each  side  of  fQ.  Figures  3.11 
and  3.12  show  log  | B(f q ,f ) j 2  for  one  sample  and  then  averaged  over 
64  samples  with  no  noise.  In  Figures  3.13  and  3.14  Gaussian  noise 
was  added  at  the  previous  levels  and  averages  over  64  samples  were 
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taken 


Log  |B(f0,f)| 


Log  IB(fo.f)  t 


3.4 


Discussion 


Comparison  of  the  noise-free  power  spectrum  (Figure  3.4), 
the  diagonal  bispectrum  B(f,f)  (Figure  3.8),  and  B(fQ,f) 

(Figure  3.12),  provides  some  measure  of  the  utility  in  signal 
classification  provided  by  the  magnitude  of  the  bispectrum.  The 
spectrum  shows  little  more  than  a  sequence  of  harmonic  lines  at 
fQ,  3fQ,  . ,.(2n  +  1)  fQ.  The  magnitude  of  the  diagonal  blspectrum, 
and  even  more  so  the  magnitude  of  B(f0,f),  shows  evidence  of  inter¬ 
actions  between  multiples  of  f_  and  the  subharmonic  f  .  This 
kind  of  information  should  clearly  permit  discrimination  of  differ¬ 
ent  signal  sources  which  have  identical  power-spectral  signatures. 
The  phase  of  the  bispectrum  should  provide  even  more  description  if 
it  does  not  depend  on  geometry  and/or  source  operating  conditions. 

All  in  all,  the  cuts  of  B(fpf2)  we  have  taken  do  not  show 
the  blspectrum  to  be  a  very  valuable  indicator  of  the  presence  of 
multiply  harmonic  signals  in  the  presence  of  Gaussian  noise. 

Several  caveats  are,  however,  in  order: 

(1)  Our  signal  had  weak  second  harmonics  and  subhar¬ 
monics.  We  examined  spectra  and  bispectra  for  another 
value  of  the  parameter  g  for  which  strong  harmonics  of 
all  orders  were  present.  The  bispectrum  was  rich  in 


detail  and  structure ,  but  no  signals  were  observed  in 
the  presence  of  noise  using  bispectra  when  they  were 
not  also  visible  in  the  power  spectra. 

(2)  We  took  only  two  cuts  through  the  two-dimensional 
frequency  space  of  the  blspectrum  B(fpf2).  One 
should  really  display  the  whole  2  dimensional  B(f^,f2) 
surface  to  exploit  the  value  of  B(fpf2>  without 

a  priori  bias  about  which  cut  has  the  essential 
Information.  The  results  of  Kim  and  Powers  (IEEE 
Trans  Plasma  Sci.,  Vol.  PS-7,  No.  2,  June  1979, 
pp.  120-131)  support  this  view.  We  have  simply  not 
had  the  time  or  computer  expertise  to  explore  this 
aspect  of  B(f ^ ,f 2) •  It  is  not  a  formidable  task. 

(3)  Our  averaging  ensembles  had  only  64  members.  Since 
Incoherent  averaging  causes  the  standard  deviation  of 
the  bispectrum  to  disappear  only  as  (degrees  of 
freedom)-*^,  we  might  have  expected  only  a  factor  of 
8  Improvement  at  best  in  our  averaging  process.  This 
clearly  argues  for  bigger  ensembles.  We  did  not 
explore  averages  of  the  bispectrum  over  an  area  in  the 
fpf2  bispectrum  space.  This  leads  to  faster  noise 
suppression  but  might  also  lead  to  loss  of  signal 
bispectrum  unless  the  bi-phase  changes  slowly  with  its 
two  frequency  arguments. 


4.0 


BISPECTRAL  LOCALIZATION 


4.1  Introduction 

There  are  many  problems  where  phase  is  irrelevant.  An 
ocean  wave  record  taken  from  10:00  to  11:00  will  not  differ  signifi¬ 
cantly  from  one  taken  between  10:01  and  11:01.  However,  the  phase 
difference  between  two  nearby  recorders  gives  significant  informa¬ 
tion  about  wave  direction.  In  order  for  this  phase  difference  to  be 
measured  with  adequate  precision,  the  two  recorders  must  be  suffi¬ 
ciently  close  to  give  coherent  records. 

This  is  the  case  where  cross-power-spectral  analysis 
provides  information  in  phase  difference  and  coherence  between  two 
records  at  one  frequency.  In  the  case  of  auto-bispectral  analysis 
we  obtain  relative  phase  information  and  bispectral  coherence 
between  two  frequencies  from  one  record.  (The  next  step  of  cross- 
bispectral  analysis  between  frequency  f ^  in  record  1  and  f2  in 
record  2  is  not  of  interest  at  the  moment.)  We  need  to  say  what  is 
meant  by  "relative  phase"  at  two  frequencies.  The  simplest 
procedure  is  to  (i)  split  the  record  into  two  by  pass-band  filtering 
at  f^  and  fj,  (ii)  heterodyne  both  to  a  single  (possibly  zero) 
frequency,  and  (ill)  obtain  the  a  cross  spectrum  between  the  two 


records.  The  blspectrum  provides  a  systematic  way  for  doing  just 
that.  Step  (iii)  is,  of  course,  in  the  spirit  of  power  spectral 
analysis,  and  so  some  will  argue  that  the  procedure  outlined  is  just 
an  application  of  ordinary  power  spectra.  We  will  take  the  view 
that  any  analysis  involving  two  frequencies  is  distinctively  differ¬ 
ent  from  the  single  frequency  power  spectrum  analysis. 

4.2  The  Bi-Phase 

Consider  two  sources  aboard  a  submarine,  separated  by  a 
distance  D,  and  emitting  signals  s'(t)  and  s'’(t),  respectively.  At 
a  receiver  R  the  received  signals  are 

s'(t)  -  E  S,<o»)ei(“t+*')  ,  -  V  -  k(u>)r' 

u> 

s"(t)  -  E  S' ' (ui)  ei*“t+*,,)  ,  t"(u)  -  v  -  k(ai)r" 

b> 

with  k  »  b)/c  ;  here  <t>'  and  are  the  relative  phases  of  the 

two  sources  and  depend  on  u>  .  For  illustration,  ambient  noise  is 
neglected  so  that  the  receiver  hears  only  the  sum  of  the  two  signals 

s(t)  -  s'  +  s”  =  E  S(u>)  eiU>t  ,  (4.1) 
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and  it  follows  that 


We  compute  B(ft,ft)  from  (4.3),  which,  on  comparison  with  (4.4), 
yields 

♦BS  “  2*'  "  *"  "  2*’  -  ♦  "  -  2k(r*  -  r”) 

with  k  -  fl/c  . 

Figure  4.1  sketches  the  locus  of  the  hyperbolae  which  have 
equal  range  difference,  r*  -  r ' ' ,  from  the  sources  s’  and  s1'. 

These  are  expressed  in  fractions  of  the  source  separation  D.  For 
the  submarine  at  a  fixed  range  r,  changing  orientation  is  equivalent 
to  moving  around  a  circle  on  Figure  4.1.  Suppose  a  trailed  submar¬ 
ine  turns  on  its  track  (which  is  standard  operating  procedure). 

O  O 

Then  6  goes  from  0  to  180  ,  and  changes  by  4  kD.  This 

ought  to  give  an  early  warning;  the  measured  phase  shift  A$ 

0b 

yields  an  estimate  of  D  which  is  of  diagnostic  value.  A  further 
diagnostic  is  the  value  $Bg  *  2$'  -  at  mid-maneuver.  The 

time-history  of  $  _(t)  for  small  r  is  slightly  different  than  for 
large  r  which  provides  some  information  about  range.  Similar 
Information  might  be  obtained  from  an  increase  in  Doppler  as  the 
target  turns  toward  the  trailer,  assuming  that  <■>'  and  <■)' '  do  not 
change.  But  the  saae  (erroneous)  conclusion  would  have  been  reached 
from  an  increase  in  source  frequency.  However,  differential  Doppler 
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Figure  4.1  D  is  the  distance  between  the  sources  s'  and  s".  The  differential  distance  Sr  = 
r'  -  r"  is  determined  by  the  bi-phase.  The  figure  shows  lines  of  constant  Sr  which 
are  hyperbolae  x2/a2-  y2^2  =  1  with  2a  =  Sr  and  2b  =  [D2  -  (Sr)2]1'2.  For 
r',  r"»  D  the  asymptote  Sr  =  d  cos0  obtains,  showing  that  bi-phase  primarily 
determines  the  orientation  of  the  source  of  the  source-receiver  path 


measurements  on  oj*  and  o>"  ,  separately,  would  provide  the  same 
localization  information  by  making  use  of  the  change  in  inter¬ 
frequency  phase. 

4.4  Interpretation  of  Sources 

The  example  of  an  ft  generator  at  one  point  well  separated 

from  a  2ft  generator  at  a  second  point  is,  of  course,  naive.  If 

the  ft  and  2ft  frequencies  are  generated  at  both  points  but  with 

unequal  intensity,  then  there  is  still  this  kind  of  information 

contained  in  $  ,  but  the  interpretation  is  more  difficult. 

BS 

Physically  we  might  expect  that  the  submarine  is  set  into 
normal  modes,  but  that  the  spatial  distribution  for  ft  modes  differ 
from  that  for  2ft  modes.  The  interpretation  is  then  similar,  with 
D  representing  the  distance  between  the  ft  mode  centroid  and  the 
2ft  mode  centroid . 

A  fascinating  speculation  is  whether  similar  estimations 


can  be  made  from  the  broad-band  acoustic  spectrum.  This  Is  not 
Impossible.  Suppose  the  flow  past  the  bow  results  In  some  resonant 
oscillation,  and  that  these  develop  harmonics  downstream  associated 


with  hydrodynamic  non-linearities.  Then  again  we  may  have  a 
situation  such  as  the  one  modeled  by  the  simple  sources  above. 

4.5  Other  Frequency  Ratios 


One  does  not  want  to  be  restricted  to  2:1  frequency  ratios 
Thus,  one  could  go  to  the  trispectrum  for  3:1  ratios,  etc.,  but  the 
formalism  becomes  awkward.  Futhermore,  how  would  one  deal  with  a 


5:3  gear  ratio? 


BIBLIOGRAPHY  ON  THE  BISPECTRUM 


Akaike,  H.  1966.  Note  on  Higher  Order  Spectra;  Annals  of  the 
Inst,  of  Stat.  Mathe.,  18:123-26. 

Armstrong,  J.W.  1977.  Bispectral  analysis  of  meter  wavelength 
Interplanetary  scintillation;  Astron.  Astrophys.,  61:313-320. 

Aubry,  M.P.  1967.  Application  de  1 'analyse  bispectrale  a 
1' etude  de  la  diffraction,  Deuxieme  partie;  Annales  d' Astrophys lque, 
30:101-110. 

Barnett,  T.P.,  L.C.  Johnson,  P.  Naitoh,  N.  Hicks  and  C.  Nute. 
1971.  Bispectrum  analysis  of  electroencephalogram  signals  during 
waking  and  sleeping;  Science,  172:401-402. 

Bartlett,  M.S.  1967.  Some  remarks  on  the  analysis  of  time- 
series;  Biometrika,  54:25-38. 

Bor resen,  R.  1978.  Experimental  determination  of  the  quadratic 
transfer  function  governing  slowly  oscillating  phenomena  in 
Irregular  waves;  Offshore  Tech.  Conf.  Proc. ,  1:457-64. 

Brilllnger,  D.R.  1973.  An  empirical  investigation  of  the 
Chandler  Hobble  and  two  proposed  excitation  processes  (Proc.  of  the 
39th  session);  Bull.  Int'l.  Stat.  Inst.,  45:413-35. 

Brilllnger,  D.R.  1973.  A  power  spectral  estimate  which  is 
insensitive  to  transients;  Technometrics,  15:559-62. 

Brilllnger,  D.R.  1974.  An  introduction  to  polyspectra.  The 
London  School  of  Economics  and  Political  Science,  1351-74. 

Brilllnger,  D.R.  1974.  Cross-spectral  analysis  of  processes 
with  stationary  increments  including  the  stationary  G/G/*  queue; 

The  Annals  of  Probability,  2:815-27. 

Brilllnger,  D.R.  1974.  Fourier  analysis  of  stationary 
processes;  Proc.  IEEE,  62:1628-43. 

Brilllnger,  D.R.  1972.  On  the  number  of  solutions  of  systems 
of  random  equations;  The  Annals  of  Mathematical  Statistics, 
43:534-40. 


Brillinger,  D.R.  1977.  The  Identification  of  a  particular  time 
series  system.  Biometrika,  64:509-15. 

Brillinger,  D.R.  1970.  The  identification  of  polynomial 
systems  by  means  of  higher  order  spectra.  J.  Sound  Vib.,  12:301-13. 

Cartwright,  D.E.  1968.  A  unified  analysis  of  tides  and  surges 
round  North  and  East  Britain;  Phil.  Trans.  Roy.  Soc.  London, 
263:1-55. 

Chambers,  J.M.  1973.  Fitting  nonlinear  models:  numerical 
techniques;  Biometrika,  60:1-13. 

Dalzell,  J.F.  1972.  Application  of  cross-bi-spectral  analysis 
to  ship  resistance  in  waves;  Davidson  Lab,  Rpt.  SIT-DL-72-1606, 
Stevens  Inst,  of  Tech.. 

Fried,  D.L.  1979.  Angular  dependence  of  the  atmospheric 
turbulence  effect  in  speckle  interferometry;  Optica  Acta, 

26:597-613. 

Fry,  J.N.  and  M.  Seldner.  1982.  Transform  analysis  of  the 
high-resolution  shane-wirtanen  catalog:  the  power  spectrum  and  the 
blspectrum;  Astro.  J.,  259:474-481. 

Gabrielli,  C. ,  M.  Keddam  and  L.  Ralllon.  1979.  Randon 
signals:  third-order  correlation  measurement;  J.  Phys.  E:  Sci. 
Instrum.,  12:632-36. 

Garrett,  J.F.  1970.  Field  observations  of  frequency  domain 
statistics  and  nonlinear  effects  in  wind-generated  ocean  waves; 

Inst,  of  Oceanography,  Vancouver,  B.C.,  Canada. 

German,  V.Kh.,  S.P.  Levikov  and  A.S.  Tsvetsinskii.  1980. 
Bispectral  analysis  of  sea-level  variations;  Meteorologiya  i 
Gldrologlya,  (Soviet  Meteorology  and  Hydrology),  11:63-70. 

Godfrey,  M.D.  1965.  An  exploratory  study  of  the  bi-spectrum  of 
economic  time  series;  Applied  Statistlcs-J.  Roy.  Stat.  Soc.,  Series 
C,  14:48-69. 

Hasselmann,  D.E.  1978.  Hind-wave  generation  by  energy  and 
momentum  flux  to  forced  components  of  a  wave  field;  J.  Fluid  Mech. , 
85:543-72. 


Hasselmann,  K.  1966.  On  nonlinear  ship  motions  in  irregular 
waves;  J.  Ship  Research,  10:64-68. 

He Hand ,  K.N.,  K.S.  Lii  and  M.  Rosenblatt.  1977.  Bispectra  of 
atmospheric  and  wind  tunnel  turbulence;  Applications  of  Statistics, 
North-Holland  Publishing  Coapany,  223-49. 

Herring,  J.R.  1980.  Theoretical  calculations  of  turbulent 
bispectra;  J.  Fluid  Mech. ,  §7:193-204. 

Hinich,  M.J.  1979.  Estiaating  the  lag  structure  of  a  nonlinear 
tiae  series  aodel;  J.  Aa.  Statl.  Assoc.,  74:449-52. 

Hinich,  M.J.  1981.  Testing  for  gausslanity  and  linearity  of  a 
stationary  tiae  series;  VA.  Tech.,  Dept,  of  Econoaics,  Blacksburg, 
VA. 


Hinich,  M.J.  and  C.S.  Clay.  1968.  The  application  of  the 
discrete  fourler  transfora  in  the  estiaation  of  power  spectra, 
coherence,  and  bispectra  of  geophysical  data;  Rev.  Geophysics, 
6:347-363. 

Huber,  P.J.,  B.  Kleiner,  T.  Gasser  and  G.  Duaerauth.  1971. 
Statistical  methods  for  investigating  phase  relations  in  stationary 
stochastic  processes;  IEEE  Trans.  Aud.  &  Electroacou. ,  AD-19:78-86. 

Hung,  G.,  D.R.  Brlllinger  and  L.  Stark,  1979.  Interpretation 
of  kernels  II.  Saae-signed  1st-  and  2nd-degree  (aain  diagonal) 
kernels  of  the  human  pupillary  systea;  Mathematical  Biosciences, 
46:159-87. 

Kedem-Kimelfeld,  B.  1975.  Estiaating  the  lags  of  lag 
processes;  J.  Aaer.  Stat.  Assoc.,  70:603-05. 

Kim,  Y.C.,  J.M.  Beall,  E.J.  Powers  and  R.W.  Miksad.  1980. 
Bispectrum  and  nonlinear  wave  coupling;  Phys.  Fluids,  23:258-63. 

Kim,  Y.C.  and  E.J.  Powers.  1979.  Digital  blspectral  analysis 
and  its  applications  to  nonlinear  wave  interactions;  IEEE  Trans. 
Plasma  Scl.,  PS-7: 120-131 . 

Kim,  Y.C.  and  E.J.  Powers.  1978.  Digital  blspectral  analysis 
of  self-excited  fluctuation  spectra;  Phys.  Fluids,  21:1452-53. 


72 


ik*-'--'"'''  '-'l 


>:  * 


Kim,  Y.C.,  E.J.  Powers,  F.  Jones  and  R.W.  Mlksad.  1981. 

Digital  Bispectral  Analysis  of  Nonlinear  Have  Couplings  in  Fluids; 
Amer.  Soc.  of  Mech.  Eng.  publ.,  81-FE-32. 

Lil,  H-S.  1982.  Model  identification  and  estimation  of 
nongaussian  ARMA  processes;  Office  of  Naval  Res.,  Arlington,  VA, 
contract  ONR  N00014-81-K0003,  Gov't,  assessment  ADA-120662. 

Lli,  K.S.  and  K.N.  Helland.  1981.  Cross-Bispectrum  Computation 
and  Variance  Estimation;  ACM  Trans,  on  Mathe.  Software 
7:284-94. 

Lii,  K.S. ,  M.  Rosenblatt  and  C.  Van  Atta.  1976.  Bispectral 
measurements  in  turbulence;  J.  Fluid  Mech.,  77:45-62. 

Liu,  P.C.,  and  A.W.  Green.  1979.  Higher  order  wave  spectra; 
Proc.  Coastal  Eng.  Cong.,  16th  pub.  ASCE,  NY,  1:360-71. 

Lumley,  J.L.  and  K.  Takeuchi.  1976.  Application  of  central- 
limit  theorems  to  turbulence  and  higher-order  spectra;  J.  Fluid 
Mech.,  74:433-68. 

Madden,  T.  1963.  Spectral,  cross-spectral,  and  bispectral 
analysis  of  low  frequency  electromagnetic  data;  Natural 
Electromagnetic  Phenomena  Below  30  KC/S  Proc.  NATO  Advanced  Study 
Institute,  Bad  Homburg,  Germany,  429-50. 

Mark,  W.D.  1970.  Spectral  analysis  of  the  convolution  and 
filtering  of  non-stationary  stochastic  processes;  J.  Sound  Vib., 
11:19-63. 

McComas,  C.H.  and  M.G.  Briscoe.  1980.  Bispectra  of  internal 
waves;  J.  Fluid  Mech.,  97:205-13. 

Narayanan,  S.  1970.  Application  of  volterra  series  to 
intermodulation  distortion  analysis  of  transistor  feedback 
amplifiers;  IEEE  Trans.  Circuit  Theory,  CT-17: 518-27. 

Ohta,  M. ,  K.  Hatakeyama,  S.  Hiromitsu  and  S.  Yamaguchi.  1975. 

Q  unified  study  on  the  output  probability  distribution  of  arbitrary 
linear  vibratory  systems  with  arbitrary  random  excitation;  J.  Sound 
&  Vibration,  43:693-711. 


73 


Ohta,  M. ,  S.  Yamaguchi  and  H.  Iwashige.  1977.  A  statistical 
theory  for  road  traffic  noise  based  on  the  composition  of  component 
response  waves  and  its  experimental  confirmation;  J.  Sound  & 
Vibration,  52:587-601. 

Patzewitsch,  J.T.,  M.D.  Srinath  and  C.I.  Black.  1978. 

Nearfleld  performance  of  passive  correlation  processing  sonars;  J. 
Acoust.  Soc.  Am.,  64:1412-23. 

Pisarenko,  V.F.  1972.  On  the  estimation  of  spectra  by  means  of 
non-linear  functions  of  the  covariance  matrix;  Geophys.  J.  Roy  astr. 
Soc.,  28:511-31. 

Ponman,  T.  1981.  The  analysis  of  periodicities  in  irregularly 
sampled  data;  Hon.  Not.  Roy  astr.  Soc.,  196:583-96. 

Rao,  T.S.  and  M.M.  Gabr.  1980.  A  test  for  linearity  of 
stationary  time  series;  J.  Time  Series  Analysis,  1:145-158. 

Roden,  G.l.  and  D.J.  Bendiner.  1973.  Bispectra  and  cross¬ 
bispectra  of  temperature,  salinity,  sound  velocity  and  density 
fluctuations  with  depth  off  Northeastern  Japan;  J.  Phys.  Ocean., 
3:308-17. 

Rosenblatt,  M.  1978.  Energy  transfer  for  the  Burgers' 
equation;  Phys.  Fluids,  21:1694-97. 

Rosenblatt,  M.  1980.  Linear  processes  and  bispectra;  J.  Appl. 
Prob.,  17:265-70. 

Rosenblatt,  M.  and  J.W.  Van  Ness.  1965.  Estimation  of  the 
bispectrum;  Annals  of  Math.  Stats.,  36:1120-1136. 

Sasaki ,  K.  and  T.  Sato.  1979.  A  bispectral  synthesizer;  J. 
Acoust.  Soc.  Am.,  65:732-39. 

Sasaki,  K. ,  T.  Sato  and  Y.  Nakamura.  1978.  An  Effective 
Utilization  of  Spectral  Spread  in  Holographic  Passive  Imalng 
Systems;  IEEE  Trans.  Sonics  &  Ultrason. ,  SU-25: 177-84. 

Sasaki,  K. ,  T.  Sato  and  Y.  Nakamura.  1977.  Holographic  Passive 
Sonar;  IEEE  Trans.  Sonics  &  Ultrason.,  SU-24: 193-200. 

Sasaki,  K. ,  T.  Sato  and  Y.  Yamashlta.  1975.  Minimum  Bias 
Windows  for  Bispectral  Estimation;  J.  Sound  &  Vibration,  40:139-148. 


Sasaki,  0.,  T.  Sato  and  T.  Oda.  1980.  Laser  Doppler  vibration 
measuring  system  using  bispectral  analysis;  Appl.  Optics,  19:151-53. 

Sato,  T.,  T.  Kishimoto  and  K.  Sasaki.  1978.  Laser  Doppler 
particle  measuring  systme  using  nonsinusoidal  forced  vibration  and 
bispectral  analysis;  Appl.  Optics,  17:667-70. 

Sato,  T.  and  O.A.  Nawa.  1968.  Bispectrum  analyzer  using 
ultrasonic  light  deflectors;  Japan-Reports,  Int'l.  Congress  on 
Acoustics,  6:5-8. 

Sato,  T.  and  K.  Sasaki.  1977.  Bispectral  holography;  J. 

Acoust.  Soc.  Amer.,  62:404-07. 

Sato,  T.,  K.  Sasaki  and  M.  Taketani.  1980.  Bispectral  passive 
veloclmeter  of  a  moving  noisy  machine;  J.  Acoust.  Soc.  Am., 
68:1729-35. 

Sato,  T.  and  0.  Sasaki.  1978.  New  3-D  laser  doppler 
veloclmeter  using  cross-blspectral  analysis;  Applied  Optics, 
17:3890-94. 

Sato,  T. ,  K.  Sasaki  and  Y.  Nakamura.  1977.  Real-time 
bispectral  analysis  of  gear  noise  and  its  application  to  contactless 
diagnosis;  J.  Acoust.  Soc.  Am.,  62:382-87. 

Sato,  T.,  K.  Sasaki  and  M.  Nonaka.  1978.  Prototype  of 
bispectral  passive  imaging  systems  aiming  machine-system  diagnosis; 
J.  Acoust.  Soc.  Am.,  63:1611-16. 

Sato,  T. ,  K.  Sasaki  and  N.  Toyoda.  1979.  A  method  to  simulate 
noises  from  a  machine  system  under  the  progress  of  abnormality;  J. 
Acoust.  Soc.  Am.,  65:1331-33. 

Shaman,  P.  1964.  Bispectral  Analysis  of  Stationary  Time 
Series,  Scientific  Paper  18;  NY  Univ,  Bronx. 

Tick,  L.J.  1961.  The  estimation  of  "transfer  functions"  of 
quadratic  systems;  Technometrics,  3:563-67. 

VanAtta,  C.W.  1979.  Inertial  range  bispectra  in  turbulence; 
Phys.  Fluids,  22:1440-42. 

Yamakawa,  S.  1976.  Investigation  of  peculiarity  in  some  wave¬ 
forms;  Bull.  JSME,  19:29-36. 


Yamanouchi,  Y.  1970.  Application  of  the  multiple  input 
spectrum  analysis  and  the  higher  order  spectrum  to  the  analysis  of 
ship  response  in  waves;  Papers  of  the  Ship  Res.  Insts.,  33:63-103. 

Yao,  N-C.G  and  S.  Neshyba.  1976.  Bispectrum  and  Cross- 
Blspectrum  Analysis  of  Wind  and  Currents  off  the  Oregon  Coast,  I. 
Development;  Sch.  of  Oceano.,  Oregon  St.  Univ,  Corvallis,  OR. 

Yao,  N-C.,  S.  Neshyba  and  H.  Crew.  1975.  Rotary  cross¬ 
bispectra  and  energy  transfer  functions  between  non-gaussian  vector 
processes  1.  Development  and  example;  J.  Phys.  Ocean.,  5:164-72. 

Yasui,  S.  1979.  Stochastic  functional  fourier  series,  volterra 
series,  and  nonlinear  systems  analysis;  IEEE  Trans.  Auto.  Control, 
AC-24: 230-242. 


DISTRIBUTION  LIST 


Dr.  Marv  Atkins 

Deputy  Director,  Science  &  Tech. 
Defense  Nuclear  Agency 
Washington,  D.C.  20305 

National  Security  Agency 
Attn  RS:  Dr.  N.  Addison  Ball 
Ft.  George  G.  Meade,  MD  20755 

Dr.  Robert  Cooper  [2] 

Director,  DARPA 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Defense  Technical  Information  [2] 
Center 

Cameron  Station 
Alexandria,  VA  22314 

The  Honorable  Richard  DeLauer 
Under  Secretary  of  Defense  (R&E) 
Office  of  the  Secretary  of 
Defense 

The  Pentagon,  Room  3E1006 
Washington,  D.C.  20301 

Director  [2] 

National  Security  Agency 
Fort  Meade,  MD  20755 
ATTN:  Mr.  Richard  Foss,  A05 

CAPT  Craig  E.  Dorman 
Department  of  the  Navy,  OP-09 5T 
The  Pentagon,  Room  5D576 
Washington,  D.C.  20350 

CDR  Timothy  Dugan 
NF0I0  Detachment,  Suitland 
4301  Suitland  Road 
Washington,  D.C.  20390 

Dr.  Larry  Gershwin 
NI0  for  Strategic  Programs 
P.0.  Box  1925 
Washington,  D.C.  20505 


Dr.  S.  William  Gouse,  W300 
Vice  President  and  General 
Manager 

The  MITRE  Corporation 
1820  Dolley  Madison  Blvd. 

McLean,  VA  22102 

Dr.  Edward  Harper 
SSBN,  Security  Director 
0P-021T 

The  Pentagon,  Room  4D534 
Washington,  D.C.  20350 

Mr.  R.  Evan  Hineman 
Deputy  Director  for  Science 
&  Technology 
P.0.  Box  1925 
Washington,  D.C  20505 

Mr.  Ben  Hunter  [2] 

CIA/DDS&T 
P.0.  Box  1925 
Washington,  D.C.  20505 

The  MITRE  Corporation  [25] 

1820  Dolley  Madison  Blvd. 

McLean,  VA  22102 

ATTN:  JASON  Library,  W002 

Mr.  Jack  Kalish 
Deputy  Program  Manager 
The  Pentagon 
Washington,  D.C.  20301 

Mr.  John  F.  Kaufmann 
Dep.  Dir.  for  Program  Analysis 
Office  of  Energy  Research,  ER-31 
Room  F326 

U.S.  Department  of  Energy 
Washington,  D.C.  20545 


D-l 


DISTRIBUTION  LIST 
(Continued) 


Dr.  George  A.  Keyworth 
Director 

Office  of  Science  &  Tech.  Policy 
Old  Executive  Office  Building 
17th  &  Pennsylvania,  N.W. 
Washington,  D.C.  20500 

MAJ  GEN  Donald  L.  Lamberson 
Assistant  Deputy  Chief  of  Staff 
(RD&A)  HQ  USAF/RD 
Washington,  D.C.  20330 

Dr.  Donald  M.  LeVine,  W385  {3] 

The  MITRE  Corporation 
1820  Dolley  Madison  Blvd. 

McLean,  VA  22102 

Mr.  V.  Larry  Lynn 
Deputy  Director,  DARPA 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Dr.  Joseph  Mangano  [2] 

DARPA/ DEO 

9th  floor.  Directed  Energy  Office 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Mr.  John  McMahon 
Dep.  Dir.  Cen.  Intelligence 
P.0.  Box  1925 
Washington,  D.C.  20505 

Director 

National  Security  Agency 


Dr.  Julian  Nall  [2] 

P.0.  Box  1925 
Washington,  D.C.  20505 

Director 

National  Security  Agency 
Fort  Meade,  MD  20755 
ATTN:  Mr.  Edward  P.  Neuburg 
DDR-FANX  3 

Prof.  William  A.  Nierenberg 
Scripps  Institution  of 
Oceanography 

University  of  California,  S.D. 

La  Jolla,  CA  92093 

Mr.  C.  Wayne  Peale 
Office  of  Research  and 
Development 
P.0.  Box  1925 
Washington,  DC  2050b 

The  MITRE  Corporation 
Records  Resources 
Mail  Stop  W971 
McLean,  VA  22102 

Mr.  Alan  J.  Roberts 

Vice  President  &  General  Manager 

Washington  Operations 

The  MITRE  Corporation 

1820  Dolley  Madison  Boulevard 

McLean,  VA  22102 


Fort  Meade,  MD  20755 
ATTN:  William  Mehuron,  DDR 

Dr.  Marvin  Moss 
Technical  Director 
Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217 


Los  Alamo 8  Scientific  Laboratory 
ATTN:  C.  Paul  Robinson 
P.0.  Box  1000 
Los  Alamos,  NM  87545 

Mr.  Richard  Ross  [2] 

P.0.  Box  1925 
Washington,  D.C.  20505 


DISTRIBUTION  LIST 
(Concluded) 


Dr.  Phil  Selwyn 
Technical  Director 
Office  of  Naval  Technology 
800  N.  Quincy  Street 
Arlington,  VA  22217 

Dr.  Eugene  Sevin  [2] 

Defense  Nuclear  Agency 
Washington,  D.C.  20305 

Mr.  Robert  Shuckman 

P.0.  Box  8618 

Ann  Arbor,  MI  48107 

Dr.  Joel  A.  Snow  [2] 

Senior  Technical  Advisor 
Office  of  Energy  Research 
U.S.  DOE,  M.S.  E084 
Washington,  D.C.  20585 

Mr.  Alexander  J.  Tachmindji 
Senior  Vice  President  4  General 
Manager 

The  MITRE  Corporation 
P.O.  Box  208 
Bedford,  MA  01730 

Dr.  Vigdor  Teplitz 
ACDA 

320  21st  Street,  N.W. 

Room  4484 

Washington,  D.C.  20451 

Dr.  A1  Trivelpiece 
Director,  Office  of  Energy 
Research,  U.S.  DOE 
M.S.  6E084 

Washington,  D.C.  20585 

Mr.  James  P.  Wade,  Jr. 

Prin.  Dep.  Under  Secretary  of 
Defense  for  R&E 
The  Pentagon,  Room  3E1014 
Washington,  D.C.  20301 


Mr.  Leo  Young 
OUSDRE  (R&AT) 

The  Pentagon,  Room  3D1067 
Washington,  D.C.  20301 


